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ABSTRACT 

Superconducting thin-film metamaterial resonators can provide a dense microwave mode spectrum with potential 
applications in quantum information science. We report on the fabrication and low-temperature measurement of 
metamaterial transmission-line resonators patterned from A1 thin films. We also describe multiple approaches for 
numerical simulations of the microwave properties of these structures, along with comparisons with the measured 
transmission spectra. The ability to predict the mode spectrum based on the chip layout provides a path towards 
future designs integrating metamaterial resonators with superconducting qubits. 
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1. INTRODUCTION 

Metamaterials having both a negative permeability and permittivity, and thus a negative index of refraction with 
left-handed transmission properties, were first described several decades digo^ More recently there have been 
numerous investigations of a variety of counterintuitive optical properties in these systems, including cloaking 
and superlensing.l^ Research in this direction is closely connected to the field of photonic band-gap engineering.!^ 

In the microwave regime, left-handed transmission lines have been proposed and studied with one- and two- 
dimensional arrays of room-temperature lumped-element components 00 There have also been implementations 
of one-dimensional metamaterial microwave transmission lines with high-temperature superconducting films 

Since 1999, there has been substantial progress in the field of quantum coherent superconducting circuits, or 
qubits.!^ Many superconducting qubit implementations involve couplings between the qubits and superconducting 
microwave resonator circuits. In this architecture, referred to as circuit quantum electrodynamics, or cQed!^ 
the qubit behaves as an artificial atom that can couple to photons in the microwave resonant cavity, analogous 
to atomic QED with natural atoms in microwave or optical cavitiesIn cQED, the microwave resonators are 
typically either distributed coplanar waveguide cavities patterned from superconducting thin film^^ or three- 
dimensional waveguide cavities E3EI1 

A recent theoretical work suggested the possibility of implementing one-dimensional superconducting meta¬ 
material resonators in the microwave regime with high quality factors for coupling to superconducting qubits.!!^ 
Such a system could allow for the generation of large-scale entanglement between multiple photon modes in the 
metamaterial and could have potential applications in the emerging field of quantum simulation. 

Here we report on the fabrication of thin-film superconducting metamaterial resonators and test oscillators 
along with microwave measurements at millikelvin temperatures. We also present numerical simulations of the 
metamaterial circuits and compare the simulated and measured spectra. The metamaterial architectures that 
we present here are compatible with future circuit designs for exploring the coupling to superconducting qubits. 
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Figure 1. Circuit schematics for (a) right-handed and (b) left-handed lumped-element lines. aS' 2 i spectra simulated in 
AWR Microwave Office for (c) right-handed and (d) left-handed transmission-line resonators. 


2. DESIGN AND FABRICATION 

2.1 Metamaterials with microwave circuit components 

Although typical planar cQED implementations involve distributed coplanar waveguide (CPW) resonators fab¬ 
ricated from superconducting thin films on silicon or sapphire substratesa resonator could also be formed 
from a lumped-element transmission line, consisting of a one-dimensional chain of unit cells each with a series 
inductor Li and a capacitor Ci to ground [Fig. [^a)].^ Such a resonant circuit would exhibit a fundamental res¬ 
onance at /o = l/2N^LiCi^ where N is the number of unit cells between the input/output coupling capacitors. 
Beyond the fundamental resonance, there are evenly spaced harmonics at integer multiples of /o, as shown in a 
circuit simulation performed in AWR Microwave Office [Fig. [^c)]. The dispersion relation for such a circuit is 
an increasing function of the wavenumber, characteristic of right-handed transmission. 

As described in Refs [5, 16] PEUby interchanging the positions of the capacitors and inductors in the lumped- 
element right-handed transmission line resonator described above [Fig.j^b)], one obtains a circuit with dramati¬ 
cally different transmission properties. There is now a low-frequency infrared cut-off at fiR = 1/Air^/LiCi^ below 
which there is no transmission through the structure. For frequencies just above fjR, there is a dense forest of 
resonances that get further apart for higher frequencies. The density of resonances for / > fjR increases as more 
unit cells are added to the transmission line. As with a conventional right-handed lumped-element transmission 
line, the impedance is still given by Zq = y/Li/Ci. 

2.2 Metamaterial design and fabrication 

We design our metamaterial resonators with a target range of fjR ^ 5 GHz in order to place the densest 
portion of the metamaterial spectrum in the range of tunability for a typical superconducting transmon qubit 
The use of an asymmetric transmon with two different sizes of junction^ would allow the IR cut-off for the 
metamaterial resonator to be located between the upper and lower high-coherence flux-insensitive sweetspots. 
Our target impedance for the metamaterial lines is Zn = 50 which, when combined with fjR = 5 GHz, results 
in Cl = 300 fF and U = 0.8 nH. 

We fabricate our metamaterial structures on a Si wafer with conventional photolithographic patterning using 
a DUV stepper followed by a lift-off process with electron-beam deposition of A1 thin films, a common material 
in superconducting qubit circuits. The A1 films are 90 nm thick. For our initial metamaterial structures we have 
used interdigitated capacitors, which allow for a straightforward single-layer fabrication process, although the 
capacitors occupy a relatively large area and their physical size limits the number of unit cells that we are able 
to fit on a chip. Our present capacitor design has fingers that are 4 fim wide, with a separation of 1 fim between 
adjacent fingers and an overiap iength of 50/im [Fig. |^a, c)]. For the inductors, we use a meander-iine design 
































































Figure 2. Optical micrographs of type-A superconducting metamaterial transmission-line resonator: (a) close-up of first 
5 unit cells, (b) coupling capacitor at one end of metamaterial, (c) close-up of interdigitated fingers in capacitor, (d) 
zoomed-out image of entire chip, (e) close-up image of two inductors. 


with a linewidth of 2 /im, gaps between adjacent lines of 3 /xm, with a length for each turn of the meander line 
of 60 /rm [Fig. |^e)]. 

In order to explore the impact of the lumped-element parameter values on the resonance spectrum, we have 
fabricated two different versions of metamaterial resonators: type A and type B. In both cases, we have attempted 
to maintain the ratio of Li and Ci for Zq ^ 50 ^2 while varying their product to change fiR. On the type A (B) 
structures, each capacitor has 32 (26) pairs of fingers and each inductor has 12 (7) turns of the meander line. 

We have chosen the lumped-element parameters based on a variety of techniques for estimating the resulting 
capacitance and inductance values. For the capacitors, numerical simulations based on ANSYS Q3D yield a 
capacitance per finger pair of 8.35 fF, corresponding to = 267 fF and = 217 fF. For the inductors, from 
an analysis with Sonnet, we obtain inductance estimates Lf = 0.9 nH and Lf = 0.6 nH. We note that these 
estimates do not account for kinetic inductance. However, for our relatively wide traces on the inductors and the 
short penetration depth of Al, we estimate the kinetic inductance contribution to be less than 5% of the total 
inductance for each L/. 

Both metamaterial structures consist of 42 unit cells extending across the width of a 6 mm-wide chip [Fig.j^d)] 
with gap coupling capacitors at either end of the metamaterial [Fig. [^b)]. For the type A(B) metamaterial 
resonators, the gap capacitor is 1 (5) /im wide, in order to target a coupling capacitance of 50 (28) fF. Each 
chip has an Al ground plane surrounding the metamaterial structure and CPW leads for probing the microwave 
transmission. As is typically done in many cQED implementations, the ground plane contains a lattice of holes 
to avoid the trapping of Abrikosov vortices in regions of large microwave currents that could contribute excess 
loss.l^ 


3. DEVICE MEASUREMENTS 
3.1 Measurements of test oscillators 

Before measuring the metamaterial resonators, we have designed and tested a series of lumped-element test 
oscillators fabricated with the same process and lumped-element parameters as in the metamaterials. The test 
oscillator chips each consist of four lumped-element LC oscillators, each one capacitively coupled to a CPW 
feedline. Each test oscillator has the same inductor as in a unit cell of metamaterial A or B and the capacitor 
has the same interdigitated finger parameters as in the corresponding metamaterial line, although the capacitor 
is split in two halves that are arranged in parallel on either side of the inductor. The total number of finger 
pairs in the capacitor for oscillator number 1-4 is 29, 33, 37, or 41, respectively, thus allowing us to study the 
variation in capacitance for different numbers of fingers [Eig. ga, b)]. 

We perform our measurements of the test oscillators on an adiabatic demagnetization refrigerator (ADR) at 
a temperature of ^ 50 mK. We probe the microwave transmission 5'2i through the feedline with a vector network 
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Figure 3. Optical micrographs of test oscillator chip with parameters corresponding to type-A metamaterial resonator: 
(a) zoomed out; (b) close-up of single oscillator, (c) Plot of 1/fo vs. number of finger pairs in capacitor of each test 
oscillator along with linear fit, as described in text. 




Figure 4. Measured S 2 i{f) of metamaterial resonator (a) type A, (b) type B. Each curve is normalized with a separately 
measured baseline, as described in text. Insets show S 2 i{f) for single peak near fm for each metamaterial. 


analyzer, sending microwave signals to the chip through a coaxial driveline with 53 dB of cold attenuation for 
thermalization. The transmission signal is amplified by a HEMT mounted on the 3K plate of the ADR and 
again with a room-temperature amplifier. A cryogenic mu-metal can mounted on the 3K plate surrounds the 
sample for magnetic shielding. 

The resonance of each test oscillator results in a dip in S '21 and we fit the resonance trajectory in the complex 
plane with a standard form,then extract the resonance frequency /o and quality factor Q for each test oscillator. 
The internal loss values extracted from these fits are in the range of 10“^, consistent with loss due to two-level 
systems at interfaces and surfaces in thin-film superconducting circuits at low temperaturesEHl^Sl Upon obtaining 
/o for each test oscillator, we are able to plot I //0 vs. the number of finger pairs in the capacitor of each test 
oscillator [Fig. |^c)]. The resulting linear variation should have a slope related to the product of the specific 
capacitance per finger pair and the common inductance value for each test oscillator. The slope in Fig. [^c) 
combined with the specific capacitance per finger pair of 8.35 fF from the Q3D numerics, corresponds to an 
inductance of Lf = 1.15 nH. A similar analysis for the type B metamaterial parameters results in Lf = 0.7 nH. 
These inductance values are slightly larger than our Sonnet estimates from the previous section, which could be 
due in part to the lack of kinetic inductance in the Sonnet treatment. Also, since our analysis of the test oscillator 
data only gives us the product of the specific capacitance per finger pair and the inductance, it is possible that 
our capacitance values from Q3D could be off, which would result in a shift of the inductance values we extract. 
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Figure 5. Simulated S 2 i{f) using AWR Microwave Office for extracted circuit parameters as described in the text (blue 
line) and comparisons with measured spectra (black line) for metamaterial (a) type A, (b) type B. Simulated curves are 
offset by —60 (—65) dB below the normalized measured spectra for type A (B). 

3.2 Measurements of metamaterials 

We cool our metamaterial resonators on the ADR to ^ 50 mK with the same measurement configuration as 
described for the test oscillators in the previous section. Using the vector network analyzer, we probe 5'2i(/) 
from 2 — 15 GHz through each metamaterial resonator, type A and B (Fig.|^. In order to normalize the spectra, 
we measure the baseline transmission on our ADR below 100 mK on a separate cooldown with a superconducting 
feedline in place of the metamaterial resonator. Subtracting this baseline for each chip accounts for any sources 
of loss off of the chip, such as cables, connectors, and attenuators, as well as the gain from the amplifiers. 

From the measured 5'2i(/) spectra, there is a clear stop band at low frequencies, with low transmission below 
^ 4 GHz for type A and ^ 5 GHz for type B, characteristic of an IR cut-off. Beyond this point, there are many 
resonance peaks that initially get closer together for increasing frequency, reaching a minimum peak spacing of 
120 (160) MHz around 5 (6) GHz for type A (B) before spreading out to larger separations at higher frequencies. 
The 3 dB linewidths of the peaks just above fjR are 260 (280) kHz for the type A (B) structures, as shown in 
the insets to Fig. 


4. SIMULATIONS OF CIRCUIT PROPERTIES 

4.1 Simulations of metamaterials: numerical solutions to lumped-element circuit 

For our initial analysis of the measured 5'2i(/) for the two metamaterial resonators, we compare these with 
numerical solutions to the lumped-element circuit model using AWR Microwave Office (Fig.|^. In this approach, 
we take the values for Q, Li and Cc from the Q3D numerics and test oscillator analysis described in Sec. 2.2 and 
3.1. The AWR simulations clearly indicate an IR cut-off, although this is roughly 1 GHz above the first peak in 
the measured spectra for both type A and type B. Also, the peak density in the AWR curves is highest just beyond 
the IR cut-off, whereas the measured 5'2i spectra exhibit an increasing peak density that reaches a maximum 
^1 — 2 GHz beyond the frequency at which the first peak occurs. This difference between the measurements and 
the AWR simulations is likely due to deviations from the ideal lumped-element model with the standing-wave 
patterns being influenced by the distributed nature of the transmission in the physical structures. In order to 
explore this, we need another approach to the simulations that accounts for the actual circuit layouts. 

4.2 Simulations of metamaterials: finite-element approach 

We simulate our physical metamaterial circuit layouts using the Sonnet software tool, which takes the GAD 
file for our circuit design as an input and solves Maxwell’s equations to obtain the field distributions that can 
then be used to compute simulated 5'2i(/) curves. From the comparison between the measured 5'2i(/) and the 
Sonnet simulations in Fig. we see that the Sonnet spectra exhibit a first peak that is within ^0.5 GHz of the 














































































Figure 6. Simulated aS'2i for circuit layout for metamaterial (a) type A and (b) type B using Sonnet (red line) as 
described in text; comparisons with measured spectra (black line). Simulated curves are offset by —45 (—60) dB below 
the normalized measured spectra for type A (B). 

measured spectra. Also, the Sonnet spectra correctly capture the location of the maxima in the peak density 
being somewhat beyond the frequency of the first peak, in contrast to the lumped-element simulations from the 
previous section. The remaining deviations between the measured 5'2i(/) and the Sonnet spectra could be due 
to various effects, including kinetic inductance in the A1 traces, disorder in the values of Ci and Li between the 
different unit cells, and coupling to chip modes in the Si substrate. 

5. CONCLUSIONS 

5.1 Conclusions 

We have demonstrated that superconducting metamaterial resonators fabricated from A1 thin films with inter- 
digitated capacitors and meander-line inductors can exhibit a transmission spectrum with an IR cut-off and a 
dense mode spectrum of high-Q peaks. By varying the capacitor and inductor values in the metamaterial, we 
are able to modify the spectrum in a predictable way. Numerical simulations of the metamaterial structures 
agree with the measured transmission spectra reasonably well. The ability to predict the mode spectrum for a 
metamaterial transmission-line resonator from the circuit design parameters is an important stepping stone to¬ 
wards the integration of similar structures with superconducting qubits for applications in quantum information 
science, for example, following the strategy described by Egger & Wilhelm.C^ 
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